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History 


- \ 


Combustion  instability  is  an  organized, 
oscillatory  motion  in  a  combustion 
chamber  sustained  by  combustion. 

V 

J 

Cl  caused  a  four  year  delay  in  the 
development  of  the  F-l  engine  used  in 

\ 

the  Apollo  program 

>  2000  full  scale  tests 

>  $400  million  for  propellants  alone 
(2010  prices) 

_ _ _ / 


Irreparable  damage  can  occur  in 
less  than  1  second. 


Damaged  engine  injector 
faceplate  caused  by  combustion 
instability 


"Combustion  instabilities  have  been  observed  in  almost  every  engine  development 
effort,  including  even  the  most  recent  development  programs" 

—  JANNAF  Stability  Panel  Draft  (2010) 
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Pressure  Fluctuation,  kPa 


Experimental  Results 


Unsteady  pressure  for  a 
translating  test 
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Instability  Mechanism 
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Complementary  Codes 


Exercised  Code  Options: 


Laminar  Combustion  Closure 
Second  Order  Accurate  in  Time  &  Space 


Choked  inlet  slots  have  been  ignored 


LESLIE 

GEMS 

r 

Structured 

Unstructured 

Explicit  MacCormack 

Implicit  Dual-Time 

LES 

DES 

r 

2-Step  reduced  mechanism 


LES:  7.3/7.6M 


DES:  3.6/4M 


Mass 

Flow 

BC 


DES  Turbulence  Model 


f - \ 

Modified  Wilcox  k-co 

\ _ y 


&plij  k 
Odj 
ckpUjUJ 


+  ak 


(t-t  +  Ou 


f - \ 

The  eddy  viscosity  is  reduced  by  modifying  the  turbulent  length  scale 

\ _ y 


i-T  =  (^T:  CdeS^) 
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LES  Turbulence  Model 


r 

Solve  a  transport  equation  of  the  sub-grid  kinetic  energy 

v 

J 

dpk*>  dpuiV*  _  _sgsdut  (*»^)3/2  a 

(  PUT  +  ^ 

,  afcsgs  pvrR  &r~ 

dt  dxi  e'  A  dxi 

[U*  '  J 

^  *  Piy 

r 

Eddy  viscosity  is  found  using  a  constant  model 

v 

\ 

_ 

V T  =  Cv  A 


( - \ 

Standard  gradient  diffusion  hypothesis  closures 

V _ J 


sgs 


=  pi/ t 


+  1^6,1 
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Pressure,  kPa 


Unstable  Operating  Point 


r  >1 

Mean  Pressure  -  1.5  MPa 

V _ 


r 

Mean  Pressure  -  1.7  MPa 

V _ J 
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Fluctuating  Pressure 
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V  _ 


"\ 
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f  \ 

DES  reaches  a 
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V _ ) 
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Phase  Difference 


The  phase  difference  between  the  head  end  and  the 
downstream  end  is  captured 
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Integrated  PSD  Data 


Mode 

Experiment 
f,  Hz  p'ptp,  kPa 

/,  Hz 

DES 

Pptpi  kPa 

/,  Hz 

LES 

Pptp>  kPa 

i 

1324 

387.15 

1500 

316.748 

1500 

382.547 

2 

2655 

89.29 

3050 

86.677 

3050 

119.221 

3 

3979 

46.37 

4550 

50.525 

4550 

70.507 

4 

7940 

41.97 

5700 

17.583 

6100 

16.879 

2 

564.78 

471.533 

589.154 

r  \ 

Similar  frequency  predictions,  both  high  relative  to  the 

experiment 

V  J 

r  's 

Mode  4  is  under  predicted 
for  both  codes 

V  J 

r  'n 

Higher  amplitude  for  LES  for 
modes  1-3 

^  J 
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Unsteady  Flowfield  -  High  Pressure 
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Unsteady  Flowfield  -  Low  Pressure 


DES 


CH4  Mass  Fraction 


0,0  0.1  0.?  0.3  0.4  0.5  0,6  0.7  0.9  0.9 


CH4  Mass  Fraction 


Temperature,  K  ^Bl  K  [  ^^^Bl 

400.0  700.0  1000.0  1300,0  1600.0  1900.0  £200.0  £500.0  £900.0  3100,0  3400.0  400.0  700.0  tOOO.O  1300.0  1500  0  1900  0  ££00  0  £500.0  £900.0  3100.0  3400.0 
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Heat  Release  Cycle  -  Part  I 
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Heat  Release  MW/m3 

Start 

V - > 
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/ - \ 
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Heat  Release  Cycle  -  Part  II 
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Reignition  Event 


LES 

/  \ 

Qualitative  agreement  with 

reignition  behavior 

V _ _ _ 


(  \ 

Accumulated  methane  in  the 

_ / 


DES 


Pressure,  Pa  CH4  Mass  Fraction 


shear  layer  in  both  cases 
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Marginally  Stable  Operating  Point 


r  >1 
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Fluctuating  Pressure 
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Integrated  PSD  Data 


Mode 

Experiment 
/,  Hz  7yptp,  kPa 

/,  Hz 

DES 

Pptp’  kPa 

/,  Hz 

LES 

PptP>  kPa 

1 

1379 

121.17 

1600 

93.527 

1500 

141.330 

2 

3881 

5.86 

3250 

23.726 

3000 

43.698 

3 

6475 

16.03 

4050 

13.573 

4200 

9.623 

E 

143.06 

130.826 

194.651 

/  \ 

LES  Over  predicts  total 
amplitude 

\ _ _ _ ) 


Both  cases  severely  over 
predict  the  second 
mode  amplitude 


r 


Frequency  differences 

v _ / 
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Unsteady  Flowfield  -  High  Pressure 


DES 


CH„  Mass  Fraction 


0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.9  0.9 
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Unsteady  Flowfield  -  Low  Pressure 
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Heat  Release  Cycle  -  Part  I 
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Heat  Release  Cycle  -  Part  II 
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Summary  &  Conclusions 


Both  LES  and  DES  area  capable  of  simulating  self-excited  combustion  instability 


Agreement  between  the  simulations  and  experiments  for  the  unstable  case  was  good 


Cyclic  Heat  release  was 
captured 


LES  had  a  delayed 
reignition,  likely 
responsible  for  the 
higher  amplitudes 


Cyclic  Heat  release  was 
captured  along  with 
reignition  event 
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